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Abstract 

We calculate the next-to-leading order inclusive total cross sections for the associated production 
processes pp — > UXk + ^ in the Minimal Super symmetric Standard Model at the CERN LHC. 
Our results provide the theoretical predictions for the total cross sections for the above processes. 
The NLO QCD corrections in general enhance the leading order total cross sections significantly, 
and vastly reduce the dependence of the total cross sections on the renormalization/factorization 
scale, which leads to increased confidence in predictions based on these results. 
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The CERN Large Hadron Collider (LHC), with = 14 TeV and a luminosity of 100 
fb~^ per year |l^, will be a wonderful machine for discovering new physics. In so many new 
physical models, the Minimal Supersymmetric Standard Model (MSSM) is one of the 
most attractive models for the theorists and the high energy experimenters, and searching 
for supersymmetric (SUSY) particles, as a direct experimental evidence, is one of the prime 
objectives at the LHC. Therefore, a good understanding of theoretical predictions of the 
cross sections for the production of the SUSY particles is important. The cross sections for 
the production of squarks and gluinos were calculated at the Born level already many years 
ago 0. To date, the productions of gluinos and squarks 1^, top squarks 0], sleptons 
P, 0] and gauginos [7] at the hadron colliders in the next-to-leading order (NLO) also have 
been studied. And recently, the NLO SUSY-QCD analysis of the associated production of a 
gaugino (x) with a gluino (g) at the Tevatron and the LHC has been presented in Ref. j^. 

In this Letter we report the NLO QCD (including SUSY QCD) calculation of the associ- 
ated production of top squarks (stops) and charginos at the LHC. Similar to pp ^ gb ^ tH~ 
jioj l , which is expected to be a dominate process for the charged Higgs boson production at 
the LHC, the associated production pp ^ gb ^ UXk niay be also the dominate process for 
single top squark or chargino production at the LHC. This is due to the following reasons: 
first, the large top quark mass in stop mass matrix can lead to strong mixing, and induce 
large mass difference between the lighter mass eigenstate and the heavier one, which means 
that the phase space for the lighter stop will be great and benefit its production; second, 
besides containing a strong QCD coupling between the incoming partons, this process also 
includes an enhanced effect from the Yukawa coupling in the vertex 6 — — of the final 
states. For simplicity, in our calculation, we neglect the bottom quark mass except in the 
Yukawa coupling. Such approximations are valid in all diagrams, in which the bottom quark 



appears as an mitia 
(ACOT) scheme 



state parton, according to the simplified Aivazis-Collins-Olness-Tung 
llj . However, it was pointed out in Ref. jl^ that the approximations of 
the hard process kinematics and the introduction of conventional bottom quark densities 
will give rise to sizable bottom quark mass and kinematical phase space effects, and may 
overestimate the inclusive cross 

tom parton approach is still valid if we choose the factorization scale below the average final 
state mass: /i/ ~ Crriav = C{mi^ +m--)/2 with C ~ (1/4, 1/3). Thus, in this Letter, we 
choose /i/ = mav/3 when we use the bottom parton approximations. 
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The leading order (LO) associated production of stops and charginos proceeds through the 
subprocess g{pa)b{ph) U{pi)Xki.'P2) with an s-channel and a t-channel. The LO squared 
matrix element in n = 4 — 2e dimensions, which has been summed the colors and spins of 
the outgoing particles, and averaged over the colors and spins of the incoming ones, is given 
by 



with 



s = iPa+ Pbf, sa = s- ml - m- 

* A. 



k 



t= (Pa- ti=t-ml, t2=t- ml 

* A. 



fc 



U= {Pa-P2) , Ui=U- m- , M2 = M - • (2) 

* ^k 

In Eq. (HI), and k\f^ are the left- and right-handed coupling constants of the vertex b—ii~-XkJ 
respectively, and are defined as follows: 

^ifc = -gRliVki + fK ^ — ^- — ^-Ri2^fc2; = fK ^ — ^ — -RliUk2- (3) 

^yZmw sm p y/2mw cos p 

Here the angle /? is defined by tan p = V2/V1, the ratio of the vacuum expectation values of the 
two Higgs doublets. Matrices U/V and are the chargino and top squark transformation 
matrices from interaction to mass eigenstates jl^ , respectively. 

The NLO corrections to the cross sections can be separated into the virtual corrections 
arising from loop diagrams of colored particles and the real corrections arising from the 
radiations of a real gluon or a massless (anti) quark. The virtual corrections consist of the 
interference of the LO amplitude M^f with the one-loop amplitudes M^^ containing the 
counterterms and the self-energy, vertex and box diagrams. We carried out the calculation 
in t'Hooft-Feynman gauge and used the dimensional regularization in n = 4 — 2e dimensions 
to regularize the ultraviolet (UV), soft infrared and collinear divergences in the virtual loop 
corrections. However, this method violates the supersymmetry. In order to restore the 
supersymmetry the simplest procedure is through finite shifts in the quark- squark- chargino 
couplings As for the Dirac matrix 75, we deal with it using the "naive" scheme, in which 
the 75-matrix anticommutes with the other 7^-matrices. This is a legitimate procedure at the 
one-loop level for anomaly-free theories jl^. The QCD coupling constant Qs is renormalized 



in the modified minimal subtraction (MS) scheme except that the divergences associated 



with the top-quark and colored SUSY particle loops are subtracted at zero momentum 



lll2 



I, |l7[ . The other renormalization constants are fixed by the on-mass-shell renormalization 
scheme 18|, and the renormalization constant of the stop mixing angle is fixed as shown 
in Ref. 19|. After renormalization, is UV-finite, but it still contains the infrared (IR) 
divergences which can be expressed as 



MlliR = + e)(^n^ + ^)Ml (4) 

zTT see 



,47r/i2 , AY. 



where 



Ar = -^, Ar = -^--ln4 + 31n4-lln^. (5) 

2 3 ' ^ 6 3 m? ml 3 ml ^ ^ 

ti ti ti 



Here the infrared divergences include the soft infrared divergences and the collinear infrared 
divergences. 

The real corrections arising from the real gluon emission gb —>■ Ux^g will produce infrared 
singularities, which can be either soft or collinear. Tly^se singularities can be conveniently 
isolated by the two cutoff phase space slicing method {2^], in which two cutoffs Sg and Sc are 
introduced, and the partonic cross section of the real gluon emission can be written as 

^il = ^i + ^fk'' + d¥- (6) 

Here the hard non-collinear part a^'^ is finite and can be numerically computed using 
standard Monte-Carlo integration techniques |21j]. The soft part contains all the soft 
infrared divergences, and is given by 

s B as r(l-e) Anfij A^ , ^1 , 4^ ^7^ 



'27rr(l-2e)' s ' 'e^ e 



with 



A', = ^, Al = ~2A^,ln6s + ^ + \ln-^-3ln-^ + Un^, 

3 6 S mf mf S mf 

ti ti 

A^, = 2AI In^ 5. - In 5. + (^2 _ 1) In 1±4 + 1 In^ <tzjl + [1 + 



'3f3 A' 7-/? 2 h ' s{^-(3)' 

ti 6 Ml 27-/? s{/3-'y) 

-2Lm '' + ^'' (8) 
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where 7 = (s + m? — m?^_)/(2s) and (3 = ^/^^ — rn^ Is. The hard colhnear part contains 
the colhnear divergences, and is given by j20| 

as r(l - e) 477 
^27rr(l-2e)^ s ' e 



dz 1 — z 

+Pgg{z,e)Gg/p{xi/z)Gb/p{x2) + {xi ^ X2)] — ( )~^dxidx2, (9) 



z z 



where Gb{g)/p{x) is the bare parton distribution function (PDF), and the unre gula ted 



splitting functions Pij{z,e) can be related to the Altarelli-Parisi splitting kernels 22] as 
Pij{z, e) = Pij{z) + ePlj{z), explicitly 

Pg,{z) = C^±^, P^^iz) = -Cf{1 - z), 

Pgg{z) = 2N[-^ + 1^ + ^(1 _ ^)], p;^(^) = (10) 
i. z z 

with Gf = 4/3 and = 3. 

In addition to the real gluon emission, other real emission corrections to the inclusive cross 
section for pp tiXk NLO involve the processes with an additional massless (anti)quark 
in the final states (g = u, rf, s, c): 

g + g^ii + xi + h, (11) 

+ + (12) 

h/h + h^ U + Xl + b/b, (13) 
q + q -^ii + Xk +b. (14) 

Since the contributions from the processes pi|) - (fT^ contain the initial state colhnear sin- 
gularities, we also need to use the two cutoff phase space slicing method [2^. However, we 
only split the phase space into two regions, because there are no soft divergences here. Thus, 
according to the approach shown in Ref. |[2Q|, the cross sections for the processes with an 



additional massless quark in the final states can be expressed as (g = u, d, s, c, b) 

do'tk'^ = ^ik{0i(^ iiXk + X)[Ga/p{Xi)Gf3/p{x2) + {Xi ^ X2)]dXidX2 

(",/3) 

+ Pga{z,e)Go,/p{xi/ z)Gh/p{x2) + {xi ^ X2)] — ( y^dxidx2, (15) 

a=q,q ^ ^ 



where 



1 + (1 _ 



= ^i^' + (1 - ^)']' W = -^(1 - (16) 

The first term in Eq. ()15|) represents the non-collinear cross sections for the four processes. 
Moreover, for tlie subprocesses gg/qq ^ tit- tiX^b {q = u, d, c, s, b), assuming m^. > ^j^^-, 
the stop momentum can approach the m^-. mass shell, which will lead to singularity arising 
from the stop propagator. Following the analysis shown in Ref . ^ \ , this problem can easily 
be solved by introducing the non-zero stop widths F^. and regularizing in this way the higher- 
order amplitudes. However, these on-shell stop contributions are already accounted for by 
the LO stop pair production with a subsequent decay into a chargino and a bottom quark, 
and thus should not be considered as a genuine higher order correction to the associated 
production of top squarks and charginos. Therefore, to avoid double counting, these pole 
contributions will be subtracted in our numerical calculations below in the same way as 
shown in Appendix B of Ref. [J]. 

After adding the renormalized virtual and real corrections, the partonic cross sections 
still contain the coUinear divergences, which can be absorbed into the redefinition of the 
PDF a. NLO. genera, called ,„a. fac.on.a.on Q. Th,s p.ocedu. 1. p.ac.ice .ea.s 
that first we convolute the partonic cross section with the bare PDF Ga/p{x), and then 
use the renormalized PDF Ga/p{x, fij) to replace Ga/p{x). In the MS convention, the scale 
dependent PDF Ga/p{x, fif) is given by [2^ 

G^/p{x,f^f) = ~'\ i^r] f -Pa,{z)G,/p{x/z). (17) 

p 6 ZTT 1 yl — Z€ j jl j J x Z 

This replacement will produce a coUinear singular counterterm, which is combined with the 
hard coUinear contributions to result in, as the definition in Ref. the 0{as) expression 
for the remaining coUinear contribution: 

da'rf = d^fk[^ Y{l- 2e) (~^)']"i^^g/p(^i' ^/^^^/p^^^' + Gy/p{xi, fXf)Gb/p{x2, fxj) 

+ [^^^-^ — + A'o^a ^ ag)]Gg/p{xi, i2f)Gb/p{x2,iJ,f) 

a=b,g ^ 

+{Xi ^ X2)}dXidX2 (18) 



where 



ATiq ^ qg) = GF{2\n5s + 3/2), (19) 



A{%g ^ gg) = 2iVln5, + {UN - 2n;)/6, 



ASC 

^0 



/3 



—Gf3/p{x/y,fif)Pa,f3{y) 



with 



y 



P'Uy)- 



(20) 
(21) 

(22) 



(23) 



Finally, the NLO total cross section for pp ^iXk factorization scheme is 



+ dxidx2[Ga/p{xi, iif)Gp/p{x2, fif) + {xi ^ X2)]afi{af3 UXk + X). (24) 

Note that the above expression contains no singularities since A2 + ^2 = and + AI + 



ATib ^ bg) 
given in Ref. 



AT{g gg) = 0. The explicit expressions of o"]^, crfl'-' and afl have been 



24|. 



We now present some typical numerical results for total cross sections for the associated 
production of top squarks and charginos at the LHC In our numerical calculations the 
Standard Model (SM) parameters were taken to be aewi^w) = 1/128, rriw = 80.419 GeV, 
= 91.1882 GeV, and nit = 174.3 GeV We use the two-loop evolution of a,(/i) Q 
{as{Mz) = 0.118), and CTEQ6M (CTEQ6L) PDFs throughout the calculations of the 
NLO (LO) cross sections. Moreover, in order to improve the perturbative calculations, we 
iake the running mass mh{Q) evaluated by the NLO formula 28| with mf,{mh) = 4.25 GeV 
3|, and make the following replacement in the tree-level couplings 

mb{Q) 



rribiQ) 



1 + Arrib '' 



(25) 



with 



Anih 



2tts nr/ \ g'^mf 

mo/i tannilmr ,mr .nia) H 7, — ^ 

377 ^'^ V fei, 62, g) i67r2m^sin2/5 



fxAtI{mi^,mi^,fx) 



- J^/iM2tan/3^[(4,)2/(m,-^,M2,/i) + i(4)2/(m^^, M^, /i)] (26) 

i=l 

where At is the soft SUSY-breaking parameter, is the Higgs mixing parameter in the 
superpotential, and 

1 



I(a, 6, c) 



,21,2 1 " , 1,2 „2i " , „2„2' ^ 



(,2_^2)(52_,2)(^^(«^ + + l°g,2^ 



(27) 



In order to avoid double counting, it is necessary to subtract these (SUSY-)QCD correc- 
tions from the renormahzation constant 6mb in the following numerical calculations. In 
the calculations of the stop decay width F^-., the two-loop leading- log relations js^ of the 
neutral Higgs boson masses and mixing angles in the MSSM were used, and the CP-odd 
Higgs boson mass m^o was fixed to 200 GeV. For the charged Higgs boson mass the tree- 
level formula was used. Other MSSM parameters are determined as follows: (i) For the 
parameters Mi, M2 and fi in the chargino and neutralino matrices, we take M2 and fi as the 
input parameters, and assuming gaugino mass unification we take Mi = (5/3) tan^ 6'vi/M2 
and rrig = [asirrig) / a2)M2 (ii) For the parameters in squark mass matrices, we assume 
Mq = Mjj = Mj^ and At = = 300 GeV to simplify the calculations. Actually, the nu- 
merical results are not very sensitive to At(^b)- Moreover, except in Fig.2, we always choose 
the renormahzation scale fir = rn^v, and the factorization scale fif is fixed to rris.v/'i- 

Fig.l shows the m^-^ dependence of the LO and NLO predictions for pp — > UXk, assuming 
/i = —200 GeV, M2 = 300 GeV and tan/5 = 30, which means that two chargino masses are 
about 182 GeV and 331 GeV, respectively, and mj^ increases from 342 GeV to 683 GeV for 
m^^ in the range 100-600 GeV. One finds that the total cross section for the ^2X2 channel 
is always smallest and less than 3 fb, but the total cross sections for other channels are 
large and range between 10 fb and several hundred fb for most values of m^^. Especially for 
the tiXi channel, the total cross section can reach 1 pb for small values of m^^ (100 GeV 
< < 160 GeV), which is almost the same as ones of top quark and charged Higgs boson 
associated production at the LHC. However, when gets larger and close to m^^, the 
total cross section for the t2Xi channel is the largest. Moreover, Fig.l also shows that the 
NLO QCD corrections enhance the LO results significantly, which are in general a few ten 
percent. The associated production of ti and Xi is the most important since the total cross 
sections are the largest for m^-^ < 400 GeV. Thus we mainly discuss this channel below. 

In Fig.2 we show the dependence of the total cross sections for the tiXi production on the 
renormalization/factorization scale, assuming fi = —200 GeV, M2 = 300 GeV, tan/3 = 30, 
m^^ = 250 GeV, and fir = f^f- One finds that the NLO QCD corrections reduce the 
dependence significantly. The cross sections vary by ±15% at LO but only by ±4% at 
NLO in the region 0.5 < fif/m^^ < 2.0. Thus the reliability of the NLO predictions has a 
substantial improvement. 

The cross sections of the tiXi production as a function of tan (3 are displayed for m^^ = 
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200, 300 and 400 GeV in Fig.3, assuming /i = -200 GeV and Ms = 300 GeV. From Fig.3, 
one can see that the cross sections become large when tan (3 gets high or low, which is due to 
the fact that for the coupling h — ti — Xi low tan (3 the top quark contribution is enhanced 
while at high tan/5 the bottom quark contribution becomes large. Fig.3 also shows that the 
NLO QCD corrections enhance the LO total cross sections, and for m^^ = 200 and 300 GeV, 
the enhancement is more significant for the medium values of tan/? than for the high and 
low ones. 

In conclusion, we have calculated the NLO inclusive total cross sections for the associated 
production processes pp tiXk ^^e MSSM at the LHC. Our calculations show that the 
total cross sections for the iiXi production for the lighter top squark masses in the region 
100 GeV < < 160 GeV can reach 1 pb in the favorable parameter space allowed by the 
current precise experiments, and besides the above case the total cross sections generally 
vary from 10 fb to several hundred fb except both m^-^ > 500 GeV and the ^2X2 production 
channel, which means that the LHC may produce abundant events of these processes, and 
it is very possible to discover these SUSY particles through the above processes in the 
future experiments, if the supersymmetry exists. Moreover, we find that the NLO QCD 
corrections in general enhance the LO total cross sections significantly, and vastly reduce 
the dependence of the total cross sections on the renormalization/factorization scale, which 
leads to increased confidence in predictions based on these results. 
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FIG. 1: Dependence of the total cross sections on m^^ for the iiXj^ productions at the LHC, 
assuming = -200 GeV, M2 = 300 GeV and tan/3 = 30. 
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FIG. 2: Dependence of the total cross sections for the tiXi production at the LHC on the renor- 
malization/factorization scale, assuming fj, = —200 GeV, M2 = 300 GeV, tan/3 = 30, m^^ = 250 
GeV, fXr = fJ-f and mav = {fnti + "^y-)/2- 
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